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ABSTRACT 

We analyse new Chandra observations of the nearest (D = 4 Mpc) Seyfert 2 active galaxy, Circinus, and 
match them to pre-existing radio, infrared and optical data to study the kpc-scale emission. The proximity of 
Circinus allows us to observe in striking detail the structure of the radio lobes, revealing for the first time edge- 
brightened emission both in X-rays and radio. After considering various other possible scenarios, we show 
that this extended emission in Circinus is most likely caused by a jet-driven outflow, which is driving shells of 
strongly shocked gas into the halo of the host galaxy. In this context, we estimate Mach numbers M ~ 2.7-3.6 
and At ~ 2.8-5.3 for the W and E shells respectively. We derive temperatures of 0.74^^ keV and 0.8 — 1.8 
keV for the W and E shells, and an expansion velocity of ~ 900-950 km s _1 . We estimate that the total energy 
(thermal and kinetic) involved in creating both shells is ~ 2 x 10 55 erg, and their age is ~ 10 6 years. Comparing 
these results with those we previously obtained for Centaurus A, NGC 3801 and Mrk 6, we show that these 
parameters scale approximately with the radio power of the parent AGN. The spatial coincidence between the 
X-ray and edge-brightened radio emission in Circinus resembles the morphology of some SNR shocks. This 
parallel has been expected for AGN, but has never been observed before. We investigate what underlying 
mechanisms both types of systems may have in common, arguing that, in Circinus, the edge-brightening in the 
shells may be accounted for by a B field enhancement caused by shock compression, but do not preclude some 
local particle acceleration. These results can be extrapolated to other low-power systems, particularly those 
with late type hosts. 

Subject headings: galaxies:active - galaxies:individual(Circinus) - galaxies:jets - shock waves - X- 
rays:galaxies - 



1. INTRODUCTION 

Recent Chandra observations of the environments of sev- 
eral radio galaxies (e.g., NGC4636, Jones et al. 2002; Hydra 
A, McNamara et al. 2000; Centaurus A, Kraft et al. 2003) 
have led to significant progress in understanding the AGN- 
driven gas outflows in these systems and the role they play in 
galaxy formation and evolution (as required by models such 
as those of Croton et al. 2006; Bower et al. 2006). We now 
know that, although the most powerful radio outflows, span- 
ning hundreds of kpc, are associated with massive elliptical 
systems, smaller structures also connected to an active nu- 
cleus can be found in a variety of systems and environments, 
including spiral and disky galaxies (e.g. Gallimore et al. 2006; 
Kharb et al. 2006; Saikia & Jamrozy 2009). 

The mechanism by which these structures are produced ap- 
pears to be similar to the one we see in the most powerful 
sources, but on a smaller scale: the relativistic plasma ejected 
by the AGN interacts with the surrounding medium, push- 
ing it out in kpc-scale radio-bright lobes, and inducing shocks 
that heat the medium to X-ray emitting temperatures (see e.g. 
Heinz et al. 1998; McNamara & Nulsen 2007; Shabala & 
Alexander 2009). The temperature structure of these bubbles 
has been analysed in a variety of systems, the most famous 
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being Centaurus A, where the bow shock has been studied 
in detail after highly detailed images of the South- West lobe 
were obtained with Chandra (see e.g. Kraft et al. 2003; Cros- 
ton et al. 2009). 

The varied morphologies of the galaxies where this mech- 
anism is found, and the fact that it is most likely episodic 
(Saikia & Jamrozy 2009; Kharb et al. 2006) make understand- 
ing the energetics involved in this process fundamental to es- 
timating its impact on galaxy evolution. The thermodynamics 
of the expanding gas provide details on the energy output from 
the AGN, the radiative timescales involved in this energy be- 
ing transferred to the ISM, and the way in which the energy 
input scales with the mass of the host galaxy, its morphology 
and the power output of the AGN, which also yields estimates 
of its lifetime. The injection rate, temperature and pressure of 
the gas determine star formation triggering and quenching in 
the central regions of the host galaxy. 

Our previous results on systems such as Cen A (Croston 
et al. 2009), NGC 6764 (Croston et al. 2008b), NGC 3801 
(Croston et al. 2007) and Markarian 6 (Mingo et al. 2011), 
show that gas outflows in these systems leave distinct sig- 
natures (thermal emission inside and in the rim of the lobes, 
synchrotron emission in very powerful shocks) whose phys- 
ical properties indicate that the radio bubbles are overpres- 
sured with respect to their surroundings (as predicted by e.g. 
Capetti et al. 1999) and are, in several cases, driving shocks 
into them. For some of these galaxies (e.g. NGC 3801, Mrk 6) 
we have been able to calculate the (substantial) fraction of the 
AGN power that is being transferred to the ISM by this mech- 
anism, setting the conditions for episodic AGN activity and 
star formation rates. Although the number of suitable sources 
for X-ray studies is limited, due to their small angular sizes, 
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Table 1 

Chandra observations of the Circinus galaxy. 



Obsid Date 


Exposure (ks) 


00356 2000-03-14 


25 


12823 2010-12-17 


160 


12824 2010-12-24 


40 



for nearby galaxies this technique is particularly fruitful, since 
with long exposures we can even resolve different regions in 
the shocks driven by the lobes, allowing us to distinguish be- 
tween thermal and non-thermal (synchrotron) emission. 

In this paper we investigate the low-power AGN in the 
Circinus galaxy. Circinus is a nearby (D — 4 Mpc) spiral 
galaxy, which exhibits a complex extended radio structure 
(Elmouttie et al. 1998a). While its nucleus and X-ray bi- 
nary populations have been studied in some depth in the past 
(Smith et al. 2000; Smith & Wilson 2001; Sambruna et al. 
2000, 2001a,b), the only preexisting X-ray observation where 
the extended emission was spatially resolved was too short to 
permit the study of the fainter emission. Our new, deep Chan- 
dra data allow us for the first time to study the X-ray coun- 
terparts of the larger-scale radio structures. Our main aim in 
the present paper is to test our new X-ray data, in conjunc- 
tion with the existing radio observations, to test the different 
scenarios that might have created the extended emission we 
observe, and to derive the physical mechanisms and energet- 
ics involved. 

Although Circinus is very close to us, the fact that it lies 
so close to our own Galaxy's plane has kept it from being 
observed in more detail in the past (rin = 5.56 x 10 21 crrT 2 , 
Dickey & Lockman 1990). Even the distance to the galaxy is 
not very well constrained, with most authors relying on the es- 
timation by Freeman et al. (1977) of 4.2^-jj Mpc, with slight 
variations. We have settled for what seems to be the most 
commonly accepted value, 4 Mpc, since this degree of uncer- 
tainty (~ 5 per cent variation in distance, ~ 10 per cent in 
luminosity) should not have an impact on our conclusions, as 
we estimate it to be below the systematics. 

2. DATA 

We observed Circinus with Chandra in December 2010, for 
a total of 200 ks. We also used a pre-existing observation from 
2000, originally analysed by Smith & Wilson (2001) (see Ta- 
ble 1), adding up a total of 225 ks. We reduced our data using 
the standard CIAO 4.2 pipeline, applying the latest calibration 
files. The observations show no trace of background flares, so 
that we could use the full exposures for all of them. 

We generated event files for each Chandra observation, 
and merged them into a single image with the CIAO tool 
merge-all. We identified the point sources with wavdetect and 
cross-checked them manually for false or omitted identifica- 
tions. We then removed the point sources and the readout 
streak (which is very prominent in the longest observation, 
running roughly in the same orientation as the galaxy's disk) 
from the merged image, to keep them from interfering in the 
analysis and to achieve the best possible characterization of 
the extended emission. Finally, we filled in the resulting gaps 
with the CIAO tool dmfilth. The resulting image can be seen 
in Fig. 1. 

We used radio maps from archival ATCA (Australia Tele- 
scope Compact Array) observations at 21, 13, and 6 cm. The 
analysis of the original data is described by Elmouttie et al. 



(1998a), who also discuss the radio properties of the source in 
detail. 

In addition to the data of Elmouttie et al. (1998a), we 
also obtained wide-band (2 GHz) ATCA data collected dur- 
ing CABB (Compact Array Broadband Backend, Wilson et al. 
2011) commissioning and as part of an ATNF Summer Stu- 
dent observation of Circinus. The CABB commissioning data 
were taken on 02 April 2009 (~8 hr, 6 cm only) using a com- 
pact 352 m configuration and 13 April 2009 (~6 hr, 3 cm and 
6 cm) using a compact 168 m configuration. The ATNF Sum- 
mer Student observations were made on 16 January 2010 (~9 
hr, 3 cm and 6 cm) using a 6 km configuration. 

The ATCA CABB data were calibrated in the MIRIAD 
package (Sault et al. 1995) and flagged using the new 
MIRIAD task mirflag. Model fitting of the source was per- 
formed using uv-components in DIFMAP (Shepherd et al. 
1994) with an additional parameter used to model the spectral 
index of each component (thus accounting for spectral varia- 
tion of the source across the 2 GHz wide band). Both phase 
and amplitude self-calibration were performed iteratively to 
improve the calibration and subsequent models. The final 
model achieves a 1 sigma sensitivity of 50 yuJy/beam. 

We also compared our data to the results on the ionized 
gas of the galaxy (Elmouttie et al. 1998b), the HI emission 
(Jones et al. 1999; Curran et al. 2008) and the optical images 
of the centre of the galaxy obtained with the HST (Wilson 
et al. 2000). 

We used XSPEC 12.5 to fit models to the X-ray spectra. All 
the model parameters thus estimated are quoted with 90 per 
cent confidence uncertainties. 

3. ANALYSIS 
3.1. Imaging 

The reduced and merged X-ray image we obtained is pre- 
sented in Fig. 1. Overlaid over the X-ray image are the 13 
cm radio contours and the regions we used to extract the sur- 
face brightness profiles shown in Figs. 3 and 4. This im- 
age shows several extended structures, with different orienta- 
tions (see also Curran et al. 2008). Our proposed geometry 
is detailed in Fig. 2, which was drawn over the combined 
X-ray/radio/optical image, to reflect the scale, reach and ori- 
entation of each component as accurately as possible. Our 
Figure is similar to the picture proposed by Elmouttie et al. 
(1998a), but there are some differences, particularly in the 
interpretation of the different structures. The AGN and cir- 
cumnuclear star-forming ring are at the centre of the galaxy, 
the ring facing towards us. The galaxy's disk extends in the 
NE-SW direction, with an orientation of ~ 60° relative to the 
plane of the sky (Jones et al. 1999; Curran et al. 2008) and 
shows some X-ray emission, though it is fainter than the other 
structures. The NW lobe is unobscured, on the near side of the 
disk, and the SE lobe is partly obscured, showing a clear dip in 
X-ray emission in the areas covered by the disk. Both lobes 
appear to be edge-brightened, though the NW lobe is much 
brighter and somewhat smaller. The orientation of the lobes 
is roughly perpendicular to that of the disk. There is a dip in 
X-ray luminosity on both lobes behind the bright edges, and a 
rise towards the centre of both, suggesting that there is further 
structure within the lobes. We extracted surface brightness 
profiles for both lobes, to quantitatively verify the correspon- 
dence between the X-ray and radio structures, as well as to 
confirm the edge-brightening in both cases. The results are 
illustrated in Figs. 3 and 4. Near the AGN, coinciding with 
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Figure 1. Merged (new+archival ACIS-S observations) 0.4-5 keV Chandra emission with overlaid ATCA 1 3 cm contours and <x = 8 Chandra pixels (4 arcsec) 
Gaussian smoothing, displaying the extent and morphology of the extended emission. Any point sources except the nucleus and the SNR have been masked 
out. The X-ray emission coincides with the edges of the radio lobes (NW and SE structures), although there is some diffuse X-ray emission coinciding with the 
galaxy's disk (extending NE-SW). This image also shows the regions we used to extract surface brightness profiles for both radio lobes (the background was 
chosen from the circular, dashed regions). The numbers indicate the first and last regions in each strip. Scale: 19.4 pc arcsec -1 . The resolution of the 13 cm radio 
map is 1 1 .8 X 11.0 arcsec. 




Figure 2. Scheme of the different components we identify in the X-ray and radio images. 
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Figure 3. Surface brightness profiles for the W lobe of Circinus, as extracted 
from the regions in Fig. 1. The regions are 30 X 9.6 arcsec in size, and have 
an inclination angle of 6.9° . The distance is measured along the regions in 
the N-S direction, starting on the northernmost edge of the first region. The 
arrows indicate the N and S edges of the shell, and the hotspot, as illustrated 
in Fig. 2. Scale: 19.4 pc arcsec" 1 
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Figure 4. Surface brightness profiles for the E lobe of Circinus, as extracted 
from the regions in Fig. 1. The regions are 23 X 7.4 arcsec in size, and have 
an inclination angle of 30°. The distance is measured along the regions in 
the N-S direction, starting on the northernmost edge of the first region. The 
arrows indicate the N and S edges of the shell, and the hotspot, as illustrated 
in Fig. 2. Scale: 19.4 pc arcsec" 1 




the base of the NW lobe, quite strong X-ray emission coin- 
cides with a region of highly ionized gas that is believed to 
be an ionization cone (Marconi et al. 1994; Elmouttie et al. 
1998b; Wilson et al. 2000; Smith & Wilson 2001). The X-ray 
emission from the galaxy's halo is very faint and not directly 
discernible from the images shown. 

The radio contours show the overlap between radio and X- 
ray emission (for a detailed analysis of the radio maps see El- 
mouttie et al. 1998a). There is faint radio emission in the NE- 
SW direction, coincident with the disk of the galaxy, while 
the radio lobes extend perpendicularly to it. The pattern of 
radio emission in the lobes matches the X-ray structures per- 
fectly: the edge-brightened emission in the lobes, the dip in 
luminosity behind the edges and the rise towards the centre 
of the lobes are spatially coincident across both wavelengths, 
except in a few points. We tested different orientations and 
distances when extracting the radial profiles in Figs. 3 and 
4. The coincidence between the radio and X-ray edge peaks 
persists across most of both lobes, though it is clearer in some 
areas. We chose regions that covered all the appropriate struc- 
tures: shells, lobe interior and hotspots, while trying to avoid 
the main emission from the ionization cone, which has no ra- 
dio counterpart, and the edges of the CCD. 

The differences in the morphology of the emission between 
the two lobes can only be partly explained by orientation ef- 
fects. The E lobe is larger, and does show X-ray emission 
outside the confines of the radio lobe edges in the regions far- 
thest from the AGN. The W lobe is smaller, more spherical, 
and its emission seems to be confined within what is seen in 
the radio maps. However, the edge of the CCD is just outside 
the edge of the W lobe in our long observations, and it is pos- 
sible that we may be missing some of the outermost edge in 
this direction, in which case both lobes could exhibit a tip of 
X-ray emission outside the radio lobe. It is possible that the 
outflow may have encountered very different environments on 
its way to its current position, which could explain the asym- 
metry between the W and E lobes. 

The W structure also shows the base of the ionization cone 
detected by Marconi et al. (1994), while no counterpart is vis- 
ible in the E, most likely due to its being obscured by the 
galaxy's disk (see Figs. 5 and 2). The Chandra images show 
three distinct tails, propagating from the AGN in the W direc- 



tion. The two brighter tails to the North have [Om] and Ha 
counterparts (see the right panel of Fig. 5 and Smith & Wil- 
son 2001 ; Wilson et al. 2000; Elmouttie et al. 1998b) while the 
southernmost edge only shows up in Ha + [Nn]. This south- 
ernmost edge could be obscured, as suggested by Elmouttie 
et al., or could be produced by a different mechanism. While 
the top tail is clearly one of the edges of the ionization cone, 
it is unclear whether the other edge is given by the central or 
the southernmost tail (Elmouttie et al. 1998b). 

3.2. Spectroscopy 

As we discussed in Section 1, the nucleus of Circinus and 
its immediate surroundings have been studied in detail in the 
past. The objective of our analysis is to study the larger scale 
X-ray emission associated with the radio lobes of Circinus. 

We extracted X-ray spectra of the radio lobes from regions 
excluding the AGN and circumnuclear emission, as well as 
the (presumably) photoionized plume in the W lobe, to avoid 
contamination. Our W region also avoids most of the area 
covered by the part of the disk that lies behind the lobe, thus 
minimising any possible contamination from it in the result- 
ing spectrum. The spectra were binned to 15 counts per bin 
to allow Gaussian statistics to be used. We also binned the 
spectra to 20 counts per bin to check the robustness of the fits. 

We used regions that cover most of the radio lobes, exclud- 
ing the circumnuclear emission and any point sources, after 
verifying that the spectrum of the gas inside the lobes does not 
differ from those of the gas in the shells, presumably due to 
the low surface brightness emission from the shells contribut- 
ing in these regions. We worked with a number of different 
source and background regions to minimize contamination. 
The spectra of the W lobe are fairly consistent regardless of 
our choice of source and background, but the E lobe was prob- 
lematic in this respect: we found that the overall shape of the 
spectrum was independent of our choice of regions, but not 
entirely consistent with any simple single or double compo- 
nent model. As discussed below, this anomalous shape may 
be caused by a higher contribution of star-related emission 
from the galaxy's disk. The statistics on the E lobe are also 
poorer than those of the W counterpart, most likely due to 
the higher obscuration, making it difficult to mask out bright 
regions and still obtain a good fit. 
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Figure 5. Left: The central regions of Circinus as seen in our merged, masked Chandra image, with <x = 2 pixel (1 arcsec) Gaussian smoothing. The image 
shows the AGN, the circumnuclear star-forming ring and the ionization cone at the base of the W radio lobe. Right: composite image showing the HST data 
on [Om] emission (red), and F814W filter (broad I) continuum (green) from Wilson et al. (2000), and the H„ emission (blue) from Elmouttie et al. (1998b), 
kindly provided by Bi-Qing For, illustrating the correspondence between the structures we observe in X-rays and those visible in the optical. The arrows point to 
matching structures: the N tail of the ionization cone (1,2), the S tail (3), which is obscured in the optical images, and parts of the circumnuclear ring (4, 5, 6). 
See Fig. 2 for more details on the geometry around these structures. Scale: 19.4 pc arcsec 

We fitted to our spectra a model consisting of local Galactic 
absorption (wabs) and emission from a collisionally ionized 
gas (apec) (see Figs. 6 and 7). We approached the analysis 
by fitting the model both to the individual spectra together 
and to the coadded spectrum (generated with the CIAO tool 
combine spectra). Only the former approach is shown in the 
Figures, but the results are consistent for both methods. Given 
that the abundance is not well constrained (especially for the 
E lobe, where we have problems with both multiphase gas 
and covering absorbers), and if left free during the fitting the 
values tend to become implausibly low, we assumed a fixed 
abundance of 0. 15, compatible with the lower limit set by the 
abundances in the IGM (see e.g. Danforth et al. 2006). Abun- 
dances higher than 0.2 do not affect the estimated values for 
the temperature, but the fit statistics are much worse. See the 
last paragraph of Section 4.2 for a discussion on the implica- 
tions of a different value of Z. 

We found that the emission in the W lobe is very well fitted 
with this single temperature model, although there are some 
residuals, especially at high energies, most likely caused by 
contamination from the AGN PSF We obtained a best-fitting 
temperature of 0.74^^ keV, and a reduced^- 2 of 1.09. See 
Fig. 6 and Table 2 for details. 

The E lobe, however, proved to be more problematic. This 
area is partially covered by the galaxy's disk, and the spec- 
trum (see Fig. 7) is not successfully modelled by a single 
temperature model (reduced^ 2 ~ 1.5). Adding a powerlaw 
or a second thermal component does improve the fit, albeit 
slightly (reduced x 1 ~ 1 -2) and most of the residuals persist. 
Adding an extra absorbing column does not improve the fit. 
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Figure 6. Spectral fit to the combined spectrum of the W lobe of Circinus. 
Model: wabs*apec, kT = 0.74+g-j*? keV. 

The poor fit could be caused by emission from an unresolved 
X-ray binary population and massive stars in the disk, or un- 
avoidable contamination from the nearby, strong AGN. 

XSPEC determined a best fitting temperature for the E lobe 
of kT = 1.6*2 ■$ keV, but it is not well constrained. A de- 
tailed analysis shows local minima when the model assumes 
kT = 0.8 and kT - 1.8 keV. Thus we used these values as a 
lower and upper limit, respectively, for the shell calculations 
displayed in Table 2. 

The disk over the E lobe is problematic not only because 
it may be introducing some contamination, it might also be 
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Figure 7. Spectral fit to the combined spectrum of the E lobe of Circinus. 
Model: wabs*apec, kT = 1.6+^ keV. 

masking an E counterpart to the W ionization cone, which 
could also be contributing to our spectra. Given the different 
morphologies of the two lobes, it is also possible that the tem- 
perature structure may be different for the E structure. We will 
discuss the nature of the emission in the lobes and its structure 
in the following Sections. 

We also carried out a fit over the regions corresponding to 
the disk (NE-SW), and found that a single temperature model 
works remarkably well, with kT — 0.69^ keV. This tem- 
perature is very similar to the one we found for the W shell, 
but not so well constrained (mainly due to a higher fraction of 
contamination from the AGN). 

To test whether there are any changes in the nature of the 
emission in the brighter areas within the lobes, we extracted 
spectra from these regions. A fit to a single powerlaw model 
gave very poor results, since the thermal emission from the 
shells is still present. We then fitted to the spectra a model 
equal to the one we used for each shell plus either a power- 
law or a second thermal component. The addition of a second 
component improves the fit; however, the poor statistics (the 
brightest individual spectrum has just ~100 counts) do not al- 
low us to distinguish between a thermal and a non-thermal 
model in this case. Although the fit is slightly better with the 
non-thermal model, the difference is not statistically signifi- 
cant. We will discuss the nature of this emission in Section 
4.3. 

4. DISCUSSION 

In this Section we discuss the possible scenarios that may 
be causing the emission we observe. In Section 4.1 we dis- 
cuss the possibility that the large-scale emission is caused 
by photoionization from the AGN. In Section 4.2 we argue 
that what we are observing is, instead, a shock. As possible 
causes behind this shock, we contemplate the possibility of a 
starburst-driven superbubble (Section 4.3) or a radio jet-ISM 
interaction (Section 4.4). 

4.1. Photoionization 

Given that the edges of the ionization cone coincide with 
the base of the W shell almost exactly (see Fig. 5 and the 
discussion at the end of Section 3.1), it is possible that some 
of the line emission studied by Marconi et al. (1994) may be 
caused by collisional ionization rather than photoionization. 
However, without [Oin], [Nn] and Ha maps with a large field 



of view and good resolution to measure the intensities and ra- 
tios, it is hard to tell. Given the intensity of the AGN and 
the light from the starburst ring, photoionization must indeed 
play a role in the emission we see in this region. However, 
it is extremely unlikely that photoionization from the central 
AGN could explain the edge-brightened, large-scale emission 
we are observing further out. We calculated the extent of 
photoionization caused by the AGN following the methods of 
Wang et al. (2009) and the models from Kallman & McCray 
(1982). Given that: 



? nR 2 



(1) 



where £ is the ionizing parameter, which we assume to be be- 
tween 10 and 100 (reasonable values for AGN outflows, as 
demonstrated by Ogle et al. 2003, see also Wang et al. 2009). 
n ~ n e is the particle density (which we assume to be between 
10~ 2 and 10~' cm 3 ), R is the distance (in cm) and L is the 
unabsorbed X-ray luminosity of the AGN required to produce 
the observed ionization factor (we assumed L = 1 .6 x 10 42 erg 
s _1 , after de Rosa et al. 2012). In the case of Circinus, this 
implies that the AGN would be able to weakly ionize (£ ~ 10) 
gas up to ~ 700 pc away, but strong ionization (£ ~ 100) 
would only be produced within the inner ~ 200 pc. Given 
the relative low luminosity of the AGN (Yang et al. 2009) and 
the dense environment surrounding it, this is consistent with 
what is observed for the ionization cone, but shows that the 
edge-brightened emission in the shells cannot be accounted 
for by photoionization. Given that the X-ray emission at the 
edges of the radio lobes extends up to 2-3 kpc in projection 
(for the W and E shell respectively), even the possibility that 
the AGN was more powerful 10 3 -10 4 years ago makes pho- 
toionization unrealistic to explain the level of emission we 
are detecting at these distances, requiring densities an order 
of magnitude smaller than those we measure in these regions 
to achieve weak photoionization. 

While there are no high resolution, large scale [Om] or Ha 
maps that cover all of the ionization cone, the existing images 
(Elmouttie et al. 1998b; Wilson et al. 2000) do show other re- 
gions of the galaxy that coincide with part of the lobes, and 
there seems to be no enhancement in emission in these re- 
gions, matching the radio and X-ray data, while this enhance- 
ment is observed in the ionization cone. From a morpholog- 
ical perspective, since the base of the W shell coincides with 
the ionization cone (Marconi et al. 1994), and given the results 
above, we can infer that there is a transition between predomi- 
nantly photoionized and shock-ionized gas at these distances, 
which could be further investigated with large scale [Oiii] and 
Ha maps and more X-ray observations. 

For consistency, we extracted a spectrum of the two north- 
ernmost tails of the ionization cone (see left panel of Fig. 
5) in our longer exposure, covering distances of 150-500 pc 
from the AGN. We fitted the resulting spectrum to a photemis 
model (part of WARMABS, an analytical XSPEC model im- 
plemented from XSTAR) and foreground Galactic absorption. 
We assumed a gas density of 10 4 cirT 3 and an ionizing pow- 
erlaw with r ~ 1.6 (de Rosa et al. 2012). We found that 
photoionized emission alone does not fully account for the 
spectral shape, yielding a reduced^ 2 ~ 2.5. This may be due 
to the less than ideal background subtraction (a local back- 
ground region that can compensate for the AGN contamina- 
tion, without including the readout streak or emission from 
the star-forming ring will still include thermal emission from 
the galaxy's disk) or to a more complex spectral structure de- 
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rived from the large volume covered by our region (a require- 
ment to obtain the necessary statistics for the fit). By testing 
the addition of different continua to improve the fit, we ob- 
served that the ionization parameter remained stable, at values 
of log(^) ~ 1.7-2.0, which supports our assumptions in the 
previous paragraph. We also found very low abundance con- 
straints (~ 10 -2 ) on O, Mg and Si, and higher values (~ 0.1- 
0.3) for Fe and Ne. We were unable to constrain the abun- 
dances for the other elements, but increasing them over our 
initial assumption (0.150) increases the^ 2 . 

We also studied the likelihood of detecting an identical 
structure in the E lobe, behind the disk of the galaxy. We 
estimated the countrate we would expect in our image from 
this structure with the PIMMS tool, assuming it to have the 
same unabsorbed flux as our W ionization cone spectrum 
0^0.3- wkev,unabs = 1.5 X 10~ 12 erg cm -2 s -1 ) and a foreground 
absorbing column from the disk comparable to that of the Cen 
A dust lane (N H ~ 1.5 X 10 22 cm" 2 ). We found that we would 
only obtain ~300 counts from this structure in our image, be- 
low the emission level from the disk of the galaxy itself, and 
thus we would not be able to detect it. Although the obscura- 
tion from the disk might be lower, it is also possible that the 
E ionization cone, if present, may be fainter than its W coun- 
terpart, if the environment through which the radiation must 
pass is denser. The possibility of an ionization cone in the E 
lobe, therefore, cannot be discarded. 

4.2. A shock model 

The match between the X-ray and radio emission is striking 
in Circinus (see Figs. 3 and 4). Unlike in other, more distant 
systems, where we lack the spatial resolution to resolve these 
structures in such detail, the proximity of Circinus allows us to 
observe that the emission from the shells coincides extremely 
well in radio and X-rays. In all the radio frequencies there is 
a dip in luminosity behind the shell, and a rise in emission to- 
wards the centre of the radio lobes. Given the edge-brightened 
nature of the emission, the abrupt jump in luminosity and the 
spectral properties we derived for the gas, we consider that the 
most likely scenario is that we are observing shells of shocked 
gas around the radio lobes of Circinus. We will discuss below 
whether this shock may be starburst or jet-driven. The dip and 
rise in luminosity behind the shells suggests further complex- 
ity in their structure, since we would expect the lobes to be 
uniformly filled with gas. We will discuss the nature of this 
emission for both the starburst and jet scenarios. 

The fact that the X-ray emission does not appear to be out- 
side the radio lobes (except at the tip of the E lobe, as dis- 
cussed in Section 3.1) suggests that the same shock that is 
causing the X-ray emission in the shells (which, as detailed 
in Section 3.2, is thermal in origin) is also giving rise to en- 
hanced synchrotron emission at radio frequencies. The appar- 
ent thickness of the shells, measured from the X-ray images, 
is ~ 300 pc, and appears constant across the entire structure. 

To test the shock hypothesis we assumed two ellipsoidal 
shells (the geometry that most closely follows the shells as 
depicted in Fig. 2, with the smallest departure from spher- 
ical symmetry) with the two minor axes equal. The two- 
dimensional regions we used in our images have eccentrici- 
ties of 0.74 and 0.52 for the E and W shells respectively. The 
E structure is larger, with R inma j = 1.0 kpc, R ex , Maj = 1.4 
kpc, R in t,min = 0.7 kpc, R ex t,min = 0.9 kpc. The W shell has 

Rint,maj — 0.8 kpc, R ex t,maj — 1-1 kpc, Ri n t,min — 0.7 kpc, 

Rext,min = 0.9 kpc. These dimensions correspond to volumes 
of 7.8 x 10 64 cm 3 and 6.1 x 10 64 cm 3 for the E and W shells 



respectively. 

We derived the electron densities for both shells from the 
normalization of the thermal model (see Table 2), assuming a 
uniform filling factor. From these densities it is also possible 
to calculate the total mass of gas in the shells, the pressure, the 
total thermal energy, and the work available from the gas. We 
found the densities to be slightly higher than those we found 
for Mrk 6, but the total energy involved is lower by roughly 
two orders of magnitude, which is to be expected given the 
much lower radio power observed in Circinus. 

In order to test our hypothesis of a shock scenario we at- 
tempted to extract spectra of the gas immediately outside the 
lobes, but were unable to fit a model to the data, due to 
the poor statistics. We then generated radial profiles for the 
longest observation from concentric annuli (r min = 132 arc- 
sec, r max = 305 arcsec, Sr = 8.7 arcsec), excluding the emis- 
sion from the AGN, the lobes and the galactic disk, and a 
suitable background. The surface brightness of the gas in the 
halo is extremely low, and it barely shows up in radial profiles 
(see Fig. 8). The drop in luminosity at R~ 440 pixels (~ 4.2 
kpc) indicates the possible extent of the halo, but the errors 
are too large to consider this anything more than a qualitative 
indication. 

To better constrain the surface brightness we estimated the 
count statistics for a circular region centered in the AGN in 
the longest observation, excluding the regions coincident with 
the lobes, the disk of the galaxy, any point sources, the read- 
out streak and any CCD gaps or adjacent chips in the process. 
We used two different radii for this circular region, r = 246 
and r = 305 arcsec, to better assess the effect of the back- 
ground, which was extracted from all the remaining S3 chip. 
We detected the halo at a 3.8<x level in the 0.3-7 keV en- 
ergy range (1439 + 382 and 884 + 239 counts for the larger 
and smaller source regions respectively, corresponding in both 
cases to a surface brightness of 0.01 1 + 0.003 counts pixeF 2 ), 
and at a 5.5<x level in the 0.5-2 keV range (1259 + 245 and 
852 ±154 counts for the larger and smaller source regions re- 
spectively, corresponding in both cases to a surface brightness 
of 0.010 + 0.002 counts pixeF 2 ). Since the orientation of the 
chips in the other two observations is different, it is not possi- 
ble to apply the same regions to them, instead we extrapolated 
the statistics from the first observation to the total observed 
time. We thus obtained a l<x range for the surface brightness 
of the halo, in the 0.3-7 keV energy band, of 0.01-0.02 counts 
pixel 2 (0.04-0.08 counts arcsec -2 ). 

The constraints on the surface brightness allowed us to test 
a thermal model with different temperatures, assuming a suit- 
able extraction region outside the radio lobes (a spherical shell 
centred in the AGN, with R int - 3.68 kpc, R ext = 5.24 kpc). 
Given the range of temperatures found in the haloes of sim- 
ilar galaxies (see e.g. Tiillmann et al. 2006b; Yamasaki et al. 
2009), we derived the electron densities for kT = 0.1 keV, 0.2 
keV and 0.3 keV. 

We used our low/high values for the surface brighntess and 
the halo shell region to obtain estimated count rates for each 
temperature. We then derived the apec normalizations re- 
quired to obtain these count rates, and derived the electron 
densities and pressures from those normalizations. The results 
are displayed in Table 3. We find that the electron densities 
we derive are compatible with those estimated by Elmouttie 
et al. (1998a), based on observations of other nearby spiral 
galaxies. We also estimated the cooling time for the gas, ob- 
taining results of 0.3-3.6 x 10 8 yr for kT = 0.2-0.3 keV, and 
as low as 6.2 Myr for kT = 0.1 keV. These times are very 
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Table 2 

Results for the shells of Circinus 



Shell kT 
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"•'-(1111 


1 iq+0.24 
^ -0.13 


9 97+0.11 

z - z '-o.n 


i 07+0.09 
'•"-o.io 


4 Q8+0- 14 


4 56+0- 22 

-0.23 


3 04 +014 
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5 55+ 030 
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4.07+0'* 
8.32+0 43 


c ao+0.24 
5 - 38 -0.26 
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Note. — These results are for the temperatures discussed in the text and Z = 0.150. The columns, from left to right, show the shell, temperature, model 
normalization, electron density, mass, pressure, total thermal energy, work available from the gas filling the shells and total kinetic energy. 
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Figure 8. Radial profile of background-subtracted counts for the longest 
Chandra observation of Circinus. Scale: 1 pixel = 0.492 arcsec = 9.54 pc. 
The emission from the halo of Circinus is barely above the background level, 
but there seems to be a net count drop after R ~ 4.2 kpc (~ 440 pixels), indi- 
cating that this may be the extent of the halo. Combining the results for the 
three observations we obtained a surface brightness estimate of 4-8 X 10~ 2 
counts/arcsec 2 . 

short, since line cooling is very important at these low tem- 
peratures, suggesting that the hot gas into which the lobes are 
expanding must have a comparatively recent origin, possibly 
in a starburst and/or supernova wind. 

To verify our assumptions on the properties of the gas in 
the halo we used the relations of Tullmann et al. (2006a), 
which give the correlation between the X-ray luminosity of 
the halo and other parameters for a sample of star-forming spi- 
ral galaxies. Some of their correlations indicate a very strong 
dependence between the soft X-ray luminosity of the halo and 
parameters such as the FIR/Ha star formation rate, and be- 
tween the overall galactic X-ray luminosity and the mass of 
dust or HI in the galaxy. 

Given how obscured and nearby Circinus is, most previ- 
ous studies have focused just on the central star-forming ring 
or the supernova remnant, so that some of the parameters we 
could use to derive the luminosity of the halo are not known, 
but we were able to use the values from Curran et al. (2008); 
Jones et al. (1999); Elmouttie et al. (1998b); For et al. (2012) 
(L BM , ~ 2.9 x 10 43 erg s" 1 , L mhol ~ 9.1 x 10 42 erg s" 1 , 
M dust ~ 10 6 M Q , M stars ~ 10 11 M , SFR FIR ~3-8M 
yr l ) with Tullmann et al.'s correlations. This allowed us to 
constrain the soft (0.3-2 keV) luminosity of the halo of Circi- 
nus to 1.6-6.3xl0 39 erg s _1 from the correlations of Tullmann 
et al. Given that the extraction region we used in our analy- 
sis does not cover the full halo of the galaxy, we find that our 
results are consistent with what is expected from these corre- 
lations and the observations (see Table 3). The luminosity we 



Table 3 

Results for the external medium. 



kT 


Norm 


N e 


Lo.2-3keV 


P 


keV xl0~ 2 cm -5 


Xl0~ 2 cm" 3 


Xl0 38 erg s _1 


Xl0~ 12 Pa 


0.1 


2.80/5.70 


2.35/3.35 


0.71/1.45 


0.70/0.99 


0.2 


1.00/2.10 


1.40/2.03 


5.30/10.70 


0.83/1.20 


0.3 


0.33/0.66 


0.80/1.14 


5.04/10.07 


0.72/1.01 



Note. — The columns, from left to right, show the low/high values of the 
model normalization, electron density, luminosity and gas pressure for each 
assumed gas temperature. 

thus constrained is also compatible with the count rates we 
estimated from our images. 

Comparing the electron densities we derive for the shells 
and the external gas, we see that density ratios of ~ 1.1- 
3 can be derived from kT ext = 0.2 or kT ea = 0.3 keV, and 
smaller values (0.6-0.9) for kT exl = 0.1 keV. The small den- 
sity contrast makes the latter temperature less consistent with 
our shock interpretation, given that the density ratio p ou tlpsiwU 
tends to 4 in the case of a strong shock. Based on observa- 
tional results (see e.g. Hummel et al. 1991; Beck et al. 1994), 
the electron densities we derived from kT = 0.2 and kT = 0.3 
keV are compatible with those seen in the haloes of spiral 
galaxies. 

Although the fact that we do not directly detect the external 
gas is still a limitation, since it would allow us to rule out a 
slow shock into a cold, dense medium, or pressure conditions 
where a shock is not possible, the morphology of the X-ray 
emission, the observed properties, and the agreement of the 
extrapolated values for the external gas with what we would 
expect in our assumed scenario indicate that we are indeed 
observing a strong shock being driven into the ISM of Circi- 
nus. 

Using the observed parameters we can calculate the Mach 
number for the shock, as follows (Landau & Lifshitz): 

/ 4(r+i)(r,wr ot ,,) + (r-i) 
M ~ V 2T (2) 

where T is the polytropic index of the gas, assumed to be 5/3. 
We find Mach numbers M ~ 2.7-3.6 for the W shell (up to 
Al = 5.1 if the external temperature is closer to kT = 0.1 
keV) and M ~ 2.8-5.3 for the E shell (up to M = 7.6 for 
kT = 0.1 keV). The sound speed is ~ 264+^ km s -1 , and 
the velocity of the shells is 915^ km s _1 for the W shell and 
950^ 4 3Q km s~' for the E shell. The errors on these values 
reflect the temperature range of both the shell and the external 
gas; the uncertainty on the E shell is very large due to the 
range of possible shell gas temperatures, as indicated in Table 
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Figure 9. 1 3 cm ATCA radio map displaying the regions we used to calculate 
the spectral index for the different areas (hotspots, lobes and shells, the RMS 
error was estimated from the dashed circle South of the AGN). With this 
colour scale the full extent of the radio emission can be appreciated, showing 
the edge-brightening in the shocked shells, and the dip in luminosity imme- 
diately behind them. The resolution of the radio map is 11.8 X 11.0 arcsec. 

2. It is clear, in any case, that the velocity of the shells must 
be supersonic, and close to 1000 km s _1 , in agreement with 
what is observed in similar systems (e.g. Mrk 6, NGC 3801: 
Mingo et al. 201 1; Croston et al. 2007) 

If we assume that the shells have been expanding with a 
constant velocity, equal to the one we have derived for the W 
structure, we can constrain their age to ~ 10 6 years. The total 
energy (enthalpy of the cavity and kinetic energy from the 
shock) involved in creating both shells is ~ 2 x 10 55 erg, the 
equivalent of 10 4 supernova explosions with individual energy 
10 51 erg, assuming a 100 per cent efficiency. 

While the spectral fits of the ionization cone (see Section 
4.1) seem to support our initial assumption of Z = 0.150, this 
value is still poorly constrained. It is very likely that the areas 
of the galaxy covered by our regions contain a multi-phase 
gas, and a range of abundances, but these cannot be studied 
in detail without deeper observations. The gas density we de- 
rive from our models depends both on the abundance and the 
emission measure, so that, assuming that there are no abrupt 
variations in the abundance between the shells and the regions 
immediately outside them, for Z = 0.30© the electron density 
in the shells would be ~ 25% lower, while that of the exter- 
nal gas would go down by ~ 20%. The ratio between both 
densities, however, only depends weakly on Z, so even if our 
assumptions on the abundance were wrong, and the densities 
were smaller, our conclusions would stand. The Mach num- 
ber and the shock velocity would not change, but the total 
energy involved in the process would be slightly smaller. 

4.3. Starburst-driven bubble 

Given that Circinus is a spiral galaxy, with active star for- 
mation, it is possible that the bubbles that we are observ- 
ing perpendicular to the plane of the galaxy may be caused 
by a starburst or supernova-driven wind. As pointed out by 
Elmouttie et al. (1998a), however, the particle component 
within starburst-driven "bubbles" is expected to be a mixture 
of thermal and non-thermal electrons, with the thermal mate- 
rial causing depolarization. The W "bubble" in Circinus does 
not show a high degree of depolarization. In the E "bubble" 
the higher depolarization is most likely caused by the gas in 
the disk, which is obscuring part of the structure (for full de- 
tails on the radio polarization of Circinus see Elmouttie et al. 
1998a). 

As we saw in Section 4.2, the total energy involved in the 
creation of both "bubbles" in Circinus is ~ 2 x 10 55 erg. This 



value is similar to the one Matsushita et al. (2005) calculated 
for the superbubble observed in M82 (0.5-2 x 1 55 erg), one of 
the most powerful examples of superbubbles. The efficiency 
with which supernova explosions transfer energy to the ISM 
is quite low (3-20 per cent, depending on the density of the 
ISM and time lapse between subsequent explosions). Elmout- 
tie et al. (1998a) inferred a supernova rate of 0.02 sn yr~' from 
their radio observations. To check this value we used the far- 
infrared data from Genzel et al. (1998) (L FIR = 7 x 1O 9 L ) 
and the correlations from Condon (1992) to estimate a mas- 
sive star formation rate (S FRm>sm ~ 0.64 M yr _1 ) and de- 
rived a supernova rate of ~ 0.026 sn yr~' . With either of these 
values, the possibility of a supernova-driven wind is energet- 
ically viable in Circinus. However, the enhanced emission at 
the centre of the "bubbles" is not easy to explain, morpho- 
logically, with this model. Multiple generations of "bubbles", 
within the 10 6 yr age estimated for the larger structure, seem 
unlikely. 

Some consideration must be given to how the radio "bub- 
bles" are related to the structures closer to the AGN. Elmouttie 
et al. (1998b) suggest that there is an ongoing outflow of gas 
in the W ionization cone (with a total kinetic energy of 5 x 1 53 
erg), which could be powered by the AGN or a starburst wind. 
They suggest that the similar sizes of the Ha central ring 
and the base of the central radio emission indicate that these 
structures are related, and that an ongoing starburst wind is 
the cause of both the inner outflow and the enhanced radio 
emission at the centre of the "bubbles". Emission line ratios 
in the cone, however, suggest that the emission in the cone 
may be caused by AGN-driven photoionization. For Circinus, 
Elmouttie et al. also point out that the morphologies of the 
ionization cone and the radio "plume" are quite different, in- 
dicating that these structures may not have the same origin. 
As suggested by Schulz (1988) for NGC4151, it is possible 
that, as the shock that created the larger structures in Circinus 
made its way through the halo gas, it created a low-density 
region through which ionizing radiation, from the AGN or a 
starburst-wind, could escape, giving rise to the cone. It is thus 
entirely possible that there is an ongoing starburst wind in the 
central regions of Circinus, feeding material into the ioniza- 
tion cone, but that this is not powerful enough to be the cause 
behind the large-scale emission in the shells and at the centre 
of the "bubbles". 

The strongest argument against the starburst-driven sce- 
nario as an explanation for the large-scale emission, however, 
comes from the variations in the radio spectral index we ob- 
serve across the "bubbles". Variations in the spectral index are 
a sign of a change in the electron population; a lower value of 
a may indicate a region where particle acceleration is occur- 
ring. An increasing value of a indicates radiative losses and 
an aging electron population. We calculated the fluxes for 
the shells, the areas immediately behind the "bubbles", and 
the regions with enhanced radio and X-ray emission inside 
them (which we labelled as hotspots), as illustrated in Figure 
9. We convolved the 6 cm, 13 cm and 20 cm ATCA radio 
maps so that the spatial resolution was the same in all of them 
(20.5 x 20.5 arcsec), thus ensuring that we are measuring the 
emission from the same structures in all cases. To best es- 
timate the extent of the RMS noise we chose a background 
region over the galaxy's disk, where there are flux variations, 
rather than a region off-source. We then derived spectral index 
values for each structure from the 6 cm/13 cm and from the 
13 cm/20 cm flux measurements. The results are illustrated in 
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Table 4 

Spectral indices 
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Note. — The spectral indices for the different extended structures 
(hotspots, shells and lobes) of Circinus, have been derived from the ATCA 
6 cm, 13 cm and 20 cm maps. The regions we used to calculate the relative 
fluxes are illustrated in Fig. 9. 

Table 4. 

Although our results are limited by the very low flux levels 
and by the fact that at 6 cm the structures we are observing 
are slightly different, it is clear that the spectral index does 
vary between these three regions, and consistently so for both 
the Western and Eastern structures, in the sense that there is 
flattening (alpha decreases) in both hotspots and steepening 
(alpha increases) in both shells, even if the individual values 
are slightly different (this may be attributed to low spatial res- 
olution, disk obscuration of the E lobe, and/or different spec- 
tral indices in the pre-shock electrons for both shells). These 
results are also consistent with the spectral index map from 
Elmouttie et al. (1998a). As stated above, this symmetrical 
flattening of the spectral index at the centre of the "bubbles" 
is most easily explained if there is ongoing particle accelera- 
tion. The simplest scenario to explain this is that these regions 
are points where a radio jet interacts with the ISM, and thus 
the equivalents of the hotspots observed in more powerful sys- 
tems (see Section 4.4 for more details). It is natural to assume 
that this same jet is the one that provided the energy injec- 
tion that caused the shock we are observing. This is further 
supported by the morphology of the emission, which is coin- 
cident in X-rays and radio, and symmetrical for both lobes, 
roughly equidistant from the nucleus of the galaxy. 

4.4. Jet-driven shocks 

In light of these results we find that the scenario that best 
describes the data is that which is typically assumed for radio 
galaxies (see e.g. Scheuer 1974; Kaiser & Alexander 1997). 
In this scenario the jet emanating from the central AGN makes 
its way into the ISM of the host galaxy, transferring part of its 



energy to the surrounding gas at the termination points (or 
hotspots). This energy transfer heats up the ISM gas, creating 
a supersonic bow shock. The material that has flowed up the 
jet is observed in the form of radio lobes. 

Hotspots in X-rays are characterised by a non-thermal pow- 
erlaw spectrum, though they are generally very faint, given 
that synchrotron losses are very high at these energies. As we 
stated in Section 3.2, the statistics of the spectra of the hotspot 
regions are quite poor, and they do not allow us to distinguish 
between thermal and non-thermal emission. If the emission 
were non-thermal, it would be unlikely to have its origin in 
widespread inverse Compton scattering, since the contribu- 
tion of such a component for low-power radio galaxies is neg- 
ligible. We estimate the IC luminosity density at 1 keV to be 
~ 6 x 10 34 erg s _1 keV -1 for the case of Circinus (equivalent 
to 2-5 x 10 -5 counts s _1 in Chandra ACIS-S, depending on 
the IC powerlaw photon index), using the assumptions on the 
B field detailed in Section 4.5. Moreover, if there was a de- 
tectable contribution from inverse Compton it would appear 
uniformly across both lobes. We calculated the spectral index 
between 6 cm and 1 keV for the W hotspot, and found it to be 
~ 1 . This value is in the range of radio/X-ray ratios seen in 
FR I jets (see e.g. Harwood & Hardcastle 2012), suggesting 
that synchrotron emission could account for at least part of 
the X-ray emission we are observing. 

Although our images do not show an active jet (at either 
radio, optical or X-ray frequencies), we know that this does 
not rule out the possibility that there is one, from what is ob- 
served in powerful FR II radio galaxies, where the jets are 
rarely seen, yet the presence of hotspots clearly indicates their 
presence (e.g. Kharb et al. 2008). The jet's rest-frame emis- 
sivity may be very low in Circinus. Moreover, the orientation 
of the lobes suggests that the jet should be at a large angle to 
the line of sight, so that relativistic beaming may be reduc- 
ing the visibility of a jet that would intrinsically be visible. 
Doppler-dimmed jets and hotspots have been found in other 
Seyfert galaxies with low radio luminosities (see e.g. Kukula 
et al. 1999), illustrating how the interactions of the jet with 
different, complex environments can change the observed ra- 
dio properties of these sources. This might also explain why 
the X-ray emission in the shells extends further than the radio 
along the (presumed) direction of propagation of the radio jet. 
This geometry is expected from some jet propagation simula- 
tions (see e.g. Sutherland & Bicknell 2007). 

In this context, to further verify our shock scenario, we 
calculated the minimum internal pressure of the radio lobes 
under equipartition conditions from the radio data, fitting a 
broken power-law electron energy spectrum with p = 2 at 
low energies, steepening to p = 3 at the electron energy that 
gives the best fit to the data. We assumed ellipsoidal geome- 
try, no protons and y min = 10. We obtained a value of P ~ 6- 
9 x 10 Pa, consistent with the results obtained by Elmouttie 
et al. (1998a), and roughly 2 orders of magnitude lower than 
the pressures we derive from our X-ray data. This departure 
from minimum energy is often found in FR I radio galaxies 
(e.g. Morganti et al. 1988); we also found this effect in NGC 
3801, Mrk6 and Cen A. These results imply that there is some 
additional contribution to the internal pressure, caused by a 
large deviation from equipartition conditions, a high fraction 
of non-radiating particles, such as thermal or relativistic pro- 
tons originating from the interaction of the galaxy's gas with 
the jet, or a low filling factor. Of these, the second explanation 
is the most plausible one (see Hardcastle et al. 2007; Croston 
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Table 5 

Power scaling 

Source £i.4Gflz V,*" E,„,j i, es 

W Hz~' sr~' erg s _1 erg 

Circinus 2.2xl() 20 10 41 2 X 10 55 

CenA 1.2X10 23 10 43 

Mrk6 1.7xl0 23 1() 42 2.6 - 4.6 xl0 56 

NGC3801 1.2xl() 23 3 X 10 42 1.7 X 10 56 

Note. — The columns display the radio core luminosity, jet kinetic lumi- 
nosity, and total energy (thermal+kinetic) inferred from the radio lobes for 
different galaxies, illustrating power scaling. 

et al. 2008a). 

From the age of the shells and the total energy we can infer 
the kinetic luminosity of the (invisible) jet, which is ~ 10 41 
erg s _1 . Comparing these results to those we obtained for Mrk 
6 (Mingo et al. 201 1), NGC 3801 (Croston et al. 2007) or Cen 
A (Kraft et al. 2003) shows that these parameters scale down 
with the radio power of the parent AGN (see Table 5). 

The fact that the X-ray emission in the shells is thermal is 
expected, if our shock scenario is correct. As pointed out by 
Croston et al. (2009), it is likely that fast shocks are required 
to produce in-situ particle acceleration and X-ray synchrotron 
emission in the shocked shells, as observed in Cen A. The 
velocities we derive, smaller than those found in Cen A, and 
the unequivocally thermal nature of the emission we observe, 
suggest that the case of Circinus is more similar to that of 
NGC 3801 or Mrk 6. 

4.5. Implications for particle acceleration and parallels with 

SNR 

The fact that the radio and X-ray emission are spatially 
coincident in the shells of Circinus resembles the morphol- 
ogy that is observed in some SNR, and poses the question 
of how similar the underlying physical processes are. Given 
that we do not observe non-thermal emission in the X-rays 
in the shells of Circinus, but we do see edge-brightening in 
radio, we investigate whether this enhancement of the radio 
synchrotron emission in Circinus can be fully accounted for 
through a shock-driven enhancement of the B field strength, 
similar to what has been observed in ongoing cluster mergers 
(e.g. Markevitch et al. 2005; van Weeren et al. 2010). The 
configuration of the magnetic field in Circinus is not very 
clear. The argument of Elmouttie et al. (1998a), stating that 
the dominant component of the magnetic field in the shells 
of Circinus seems to be predominantly radial (perpendicular 
to the shock front, parallel to its direction of propagation) is 
based on polarization maps and is not very strong. Other sys- 
tems such as NGC 1068 (Wilson & Ulvestad 1987; Simpson 
et al. 2002) and NGC 3079 (Duric & Seaquist 1988) indicate 
that fields tangential to the shock front (perpendicular to the 
direction of propagation of the shock) are more often found in 
radio galaxies. Both radial and tangential B field configura- 
tions can be found on SNR (see e.g. Milne 1987), but only the 
component that is tangential to the shock front is compressed 
by it. Compressed magnetic fields have been found around 
the leading edges of FR I radio galaxies (Guidetti et al. 201 1, 
2012) from rotation measure analysis, suggesting that, given 
the right conditions (a component to the magnetic field tan- 
gential to the front and a strong shock) this phenomenon may 
be common to many low-power systems. 

To investigate whether this shock-driven B enhancement 



could be enough to account for the enhancement in radio 
emission, we must calculate whether the jump in emissivity 
(of a factor of at least 40) observed between the shells and 
the gas just outside them can be explained by this mechanism. 
The synchrotron emissivity of a powerlaw distribution of elec- 
trons is given by : 

(p-i) p+i 

J(v) = CN v-—B— (3) 

where C can be assumed constant, No is the normalization 
of the relativistic electron spectrum, and (p - l)/2 = a. If 
we assume that B is tangential to the shock front, applying 
the Rankine-Hugoniot conditions for a strong shock results in 
Bsheiis/Bext K Psheiis/Pext, which, substituting the values from 
Tables 2 and 3, results in B s heiis ~ 1-1-3B CX/ . In Eq. 3 No is 
proportional to the electron density, but given that the elec- 
tron density depends on the energy density, which scales with 
the temperature in the shock, we have that No, s heiislNo,ext K 
(PsheitelPextW shells I T ex t) p ~ l - Substituting the values from Ta- 
bles 2, 3 and 4 we can calculate the emissivity jump. Taking 
only the values for the W shell, which are better constrained, 
we see that a ~ 0.70-0.94 and T sheU ~ 0.69 - 0.81 keV. 
If T ext = 0.2 keV, using the densities from Tables 2 and 3, 
we have (p she iisl Pext) ~ 1.1-1.7, and (T she „ s IT ex ,) ~ 3.5- 
4.1, thus {J shells I J ext) ~ 7-68. If T ext = 0.3 keV we have 
(p shells I Pext) ~ 1.9-3.0, and (7 ' shea J7 \xi) ~ 23-2.1, thus 

{J shells I J ext) ~ 18-164. 

These results show that shock-compression alone could be 
enough to fully account for the enhancement in radio syn- 
chrotron emission if the conditions are favourable (a strong 
shock and a magnetic field tangential to the shock front). It is 
likely that, even in the best conditions, there is a component 
of the magnetic field that is perpendicular to the shock front. 
Our results do not preclude local particle acceleration — it is 
in fact needed if the conditions at the shock are not ideal — 
but we know that any particle reacceleration is not efficient 
enough to power synchrotron emission in the X-rays. 

When making the parallel between the shells of Circinus 
and what is observed in SNR, some caution must be exer- 
cised. In young SNR where B field amplification is observed, 
it is mediated by efficient particle acceleration (Volk et al. 
2005). It seems to be clear that, in most cases, field compres- 
sion alone cannot account for the acceleration of cosmic rays 
that is observed in SNR, which are in fact the main source of 
Galactic cosmic rays. The field amplification in SNR is much 
larger than the enhancement expected from shock compres- 
sion alone, given that the values of B in SNR tend to be of 
the order of mG (Vink & Laming 2003), 2-3 orders of mag- 
nitude higher than the values found in spiral galaxies, where 
typical values for the ISM of the disk are ~ 10 pG, (see e.g. 
Beck et al. 1996; Vallee 2011; Hummel et al. 1991; Beck 
et al. 1994). Synchrotron losses also become substantial at 
the values of B found in most SNR. Rayleigh-Taylor instabil- 
ities at the shock front (Guo et al. 2012) seem to be present 
in most SNR. While Rayleigh-Taylor instabilities may be be- 
hind some of the B enhancement in shock shells around radio 
galaxy lobes, strong instabilities are not expected to develop 
until the jet has switched off and the radio lobes have become 
buoyant (see e.g. Briiggen & Kaiser 2001). 

The fact that the spectral index in the hotspots of Circinus is 
flatter than that of the lobes suggests that there is spectral age- 
ing of the synchrotron emission, evolving radially outwards 
from the hotspots. In the model we discuss here the electron 
population in the shells, however, must come from electrons 
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from the external medium swept up in the shock. These elec- 
trons are already relativistic, as we know from synchrotron 
observations of the haloes of spiral galaxies (see e.g. Beck 
et al. 1996). We cannot directly measure the spectral index of 
the emission outside the shells, since it is too faint. However, 
we studied the spectral index in the disk, and found a value 
a ~ 1 . This value is similar to that found in the shells, hint- 
ing that this is the spectral index of the electron population 
outside the shells as well. This provides additional support 
to the idea that the halo electrons swept up in the shock are 
not strongly reaccelerated, since that would flatten their spec- 
tral index, and that the main mechanism behind the enhanced 
synchrotron emission in the shells we observe in the radio 
maps is indeed a B field enhancement. The thermal nature of 
the X-ray emission in the shells also supports the hypothesis 
of shock compression as the dominant mechanism behind the 
radio emission enhancement. 

We carried out equipartition calculations for the shells, and 
found that, under those conditions, the pressures we derive are 
over an order of magnitude below our measurements. This 
suggests that the shells are not in equipartition and that the 
contribution to the total energy density from the synchrotron- 
emitting electrons is not the dominant one. This is to be ex- 
pected, since the spatial coincidence between the X-ray emit- 
ting hot gas and the synchrotron-emitting plasma implies that 
they coexist in this region, and the former component is ex- 
pected to dominate the energy density. 

As mentioned above, compressed magnetic fields may be 
common to low-power radio galaxies. Given that field en- 
hancement through shock compression does not preempt par- 
ticle acceleration, it is likely that in most systems a combi- 
nation of both is responsible for the observed emission. An 
accurate measurement of the spectral index and emissivity of 
the electrons outside the lobes, and better X-ray statistics for 
both the shells and the external environment, would allow us 
to estimate the contributions from both processes in Circinus 
and in other systems. In the case of systems more similar to 
Cen A, where particle acceleration seems to be causing syn- 
chrotron emission in the shells, all the way to the X-rays, con- 
straints on the magnetic field in the shock, e.g. via rotation 
measure analysis, as in Guidetti et al. (2011), could be used 
for estimations on the contribution to the emission caused by 
magnetic field compression. 

Our results suggest that the parallel between the jet-driven 
shocks of low power FR I galaxies and those of SNR, where a 
similar morphology and some of the same physical processes 
can be expected, is more common than we previously thought. 
The fact that it has only been observed in Circinus is mainly 
due to the small distance between this galaxy and our own, 
which allows us to spatially resolve these structures, and ob- 
serve their characteristics with the level of sensitivity of our 
current instruments. This is particularly relevant for more dis- 
tant low-power radio galaxies in general, and those with late- 
type hosts in particular. It is also significant for the edge- 
brightening in the radio lobes of Mrk 6 (Kharb et al. 2006; 
Mingo et al. 2011), which could be explained through this 
mechanism. As pointed out in Section 4.2, the shock results 
suggest that we are facing a scenario similar to the one we 
found in NGC 3801, suggesting that some of our conclusions 
may apply to this object, as well as, possibly, to NGC 6764. 

5. CONCLUSIONS 

We have presented a detailed analysis of the X-ray emis- 
sion associated with the radio lobes of the Circinus galaxy. 



Our deep Chandra observation has allowed us to study both 
the morphological correspondence between radio and X-ray 
emission, and the X-ray spectral properties of the gas in 
these structures. We conclude that we are observing shells 
of shocked gas around the radio lobes of Circinus, expand- 
ing into the halo regions of the galaxy with Mach numbers 
M ~ 2.7-3.6 for the W shell and M ~ 2.8-5.3 for the E 
shell, consistent with the Rankine-Hugoniot conditions for a 
strong shock. 

We rule out the possibility that this emission is caused by 
AGN photoionization, and discard the scenario in which the 
radio structures are the result of star formation or a supernova- 
driven wind; instead we argue that they are caused by the in- 
teraction of a radio jet with the surrounding gas. The total 
energy (thermal and kinetic) involved in the creation of these 
shells is ~ 2 x 10 55 erg. We estimate their age to be ~ 10 6 
years. From these parameters we infer the kinetic luminosity 
of the jet, which is ~ 10 41 erg s _1 . When comparing these 
results with those previously obtained for Cen A, Mrk 6 and 
NGC 3801, we observe that they scale with the radio power 
of the source. 

The presence of two symmetrical structures resembling jet 
termination points or "hotspots", together with the slightly 
different X-ray spectra and flattening of radio spectral index 
in these regions, suggests that the radio jet may still be ac- 
tive, even if it is not visible. This suggests that "invisible" jets 
may be a more common occurrence among Seyfert and late 
type galaxies than previously thought, given that it is difficult 
to separate the different radio structures in galaxies with ac- 
tive star formation, when the radio bubbles are faint or their 
orientation is similar to that of the disk. This opens new ques- 
tions on the role of low radio luminosity, star-forming, active 
galaxies, which are often ignored in survey studies. 

In Circinus the X-ray and radio emission from the shocked 
shells are spatially coincident. This, added to the fact that 
the radio emission is edge-brightened, allows us to draw a 
parallel between these structures and what is observed in su- 
pernova remnants: the radio emission from the shells comes 
from shock-compressed or (re)accelerated cosmic ray elec- 
trons in the environment, but particle acceleration, if present 
at all, is not efficient enough to produce TeV electrons, so that 
only thermal emission is observed in X-rays. This correspon- 
dence has been long expected for AGN, given the similarity 
between the physical processes involved in SN and jet-driven 
shocks, but has not been unequivocally observed before. We 
believe this scenario may be common to other low-power ra- 
dio galaxies, and potentially very relevant to those with late- 
type hosts, but can only be observed where we are able to 
spatially resolve the different structures, something that has 
been achieved for Circinus due to its proximity to us. 
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